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Abstract 

In this letter we extend the nonlinear a model describing pure spin ladders with an 
arbitrary number of legs to the case of ladders containing a single static hole. A simple 
immediate application of this approach to classical ladders is worked out. 

PACS numbers: 75.50.Ee; 75.10.Jm; 75.10.Hk; 75.30.Ds 
Keywords: magnons; antiferromagnets; impurities. 

Corresponding author: A. R. Pereira; e-mail: apereira@ufv.br; Tel.: +55-31-3899-2988, 
Fax: +55-31-3899-2483. 



*E-Mail: apereira@ufv.br 

* E-Mail:Elisa.Ercolessi@bo.infn.it 

■'•E-mail: antpires@fisica.ufmg.br 



1 



The interest in low dimensional quantum antiferromagnets has been great ever since 
Haldane conjectured [T] that integer spin chains have a gap in their excitation spectrum 
while half- integer spin chains do not. More recently, spin ladders (two or more coupled 
spin chains) have also attracted much interest, mainly when the effects upon doping are 
considered. Indeed, when one manages to remove spins from the system (leaving holes 
behind) the existence of superconductivity is predicted [2] (and experimentally observed 
[3]). In this letter we would like to study the presence of static holes in spin ladders 
with an arbitrary number of legs by considering the nonlinear a model continum limit 
of the model. 

The continuum limit of pure antiferromagnetic Heisenberg spin ladders with arbi- 
trary number of legs has been derived by several authors [H El El [7] . Here we apply the 
approach of Ref.fl] for studying the case of ladders containing static holes. We start 
with the Hamiltonian for a pure ladder system with n legs of length Nao (ao is the 
lattice constant and N ^> n) defined as 

n N 

H = ^2J2 [ J ^a(j) ■ S a (j + 1) + J a , a+ lS a (j) ■ S a+1 (j)] , (1) 
a=l j=l 

where S a (j) are the spin operators located in the ath leg at the position j = 1, ...,N, 
while J a > and J a ,a+i > are the antiferromagnetic exchange couplings along the 
leg and rung respectively. The partition function of the above Hamiltonian in the spin 
coherent state path-integral representation is given by 




where r = it is the imaginary time variable, uj[£l a {j, r)] is the Berry phase factor and 
H(t) is obtained by replacing the operator S a (j) by the classical variable S£l a (j, r) in 
the Hamiltonian ([!]). To get the continuum limit, it is usual to assume that the dominant 
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contribution to the path integral comes from paths described by [TJ [8] 

( \ 1/2 
1 - ' laU s ^ ' j (3) 

The field (f>(J,r) is supposed to be slowly varying and the fluctuation field l a (j,r) is 
supposed to be small (l a (j,T)/S « 1). The constraint Q,^(J,t) = 1 implies </> 2 (j, r) = 1 
and 4>(j,r) ■ l a (j) = 0. Numerical works support the fact that the staggered spin-spin 
correlation length is much greater than the total width of the ladder [9j [10]. Then, 
assuming that 4>(j, r) depends only on the site index j along the legs, DelFAringa et 
al. [1] mapped the antiferromagnetic Heisenberg ladder system onto a (1+1) quantum 
nonlinear a model 



Z a = J [D<f] exp(ir[<^]) exp J dr J dx 



-{djf + v s {djf 



(4) 



where T[(f>] = (9 /Air) Jjf dr J dx<f> ■ (d T (f> x d x (f>) (with 6 = 2nS) for n odd and T[(j>} = 
for n even, reflecting the fact that, for half-spin systems, the excitation spectrum has a 
gap (is gapless) when n is even (odd). Besides, the nonlinear a model parameters, the 
coupling constant g and the spin wave velocity v s , are defined by 

a' 1 = S [ £ JaL-A , (5) 



i a,b,c 



1/2 

v s = S ( ] , (6) 



J2b,c L b 



where L a . is the inverse of the matrix 



, 4J a + J a ,a+l + Ja,a-l for CL = b 
L a,b = \ {<) 

L a ,b = J a ,a+l for I a - b |= 1 



with J a , a +i = J a -i,a and J lfi = J n , n +i = 0. 



Now we consider the system in the presence of static holes (spins removed from the 
ladder). In two spatial dimensions, one of the simplest way of studying this problem in 
the continuum limit is through a non-simply connected manifold [T2| [T3] . In this 
case a disk is removed from the magnetic plane, leaving a hole behind, and this hole is 
interpreted as a nonmagnetic impurity (or a spin vacancy) since there is no magnetic 
degrees of freedom insight it. This approach has good qualitative and quantitative 
agreement with numerical calculations 112 } [T3 l I14j. However, in the case of ladders 
as described by the nonlinear a model given by Eq. we cannot simply remove a 
disc from the space because the problem becomes essentially one-dimensional. There is 
no possibility of removing a part of the space without breaking the "effective" lattice. 
Then, within the above approach a "hole" must affect the exchange interactions J a and 
Ja,a±l- As a consequence, the matrix L a ^ becomes j-dependent. In such a way that 
the lattice is not broken for a single defect. It means that the parameters g and v s are 
now functions of the position along the ladder. Depending on the number of legs, there 
are more than one position to put a single vacancy which yields different results. Some 
examples are shown in Fig.(l). 

With the above considerations in mind, our approach for the nonlinear a model 
describing spin ladders with a static hole centered at j = xo (along the legs), a = k 
(along the rungs) are summarized as follows 



In order to explicitly define the new parameters space dependents gk(x — xo) and v s ^{x — 





with 




(9) 
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xq), we first rewrite the n x n matrix L a ^ as follows 



L l,i 
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If we place the vacancy at (J, a) = (xQ,k), the above matrix will be the same for all 
j 7^ xq, xq — 1 and equal to ([7]) . For j = xq — 1 the matrix will be again of rank n with a 
slightly different coefficient Lf- t k- F° r j = the matrix has zeroes along the k — th row 
and the k — th column. Thus we define a new matrix K a ^ of order (n — 1) x (n — 1), 
which has almost the same elements of the above matrix and without line k and column 
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k. Explicitly: 



L l,l 


Ll t 2 
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, (11) 



where K k _^ k _ x = L k _ x>k _i - J fc _i jfc and K k+Xyk+ i = L k+1)k+1 - J k+Xtk . This formula 
holds only in the region of the spin vacancy (| x — xq |< ao). Therefore we may assume 
that the coefficients g and v are given by: 



g k [x-x Q ) 



V s>k (x - Xq) 



for I x- x \> a , 



^k,b,c 



J a K t 



b.c 



1/2 



for I x- x \< a , 
for \ x- xq |> ao, 



(12) 



(13) 



The change in the parameters g k (x— xq) and v k;S (x— xq) in the zone of influence of the 
vacancy is associated with the discontinuous change in the number of legs in this region 
(see Fig.(l)). Of course, the field (j) must be continuous across the pure and impure 
regions. Note that, in principle, only the coupling constant and spin wave velocity are 
(locally) affected by the removed spin. As the Berry phases do not depend on these 
values, they are not very sensitive to the presence of the defect (see Eq. ©)■ It means 



that the ground state may not be very affected by the presence of the impurity. These 
results are in agreement with recent numerical calculations[15j, which give evidences that 
the ground state configuration of the entire ladder system is not changed significantly 
by the impurity except for the local extraction of the missing bonds. Indeed, the energy 
cost to remove a spin from a two-leg ladder with spin-1/2 and J = J is E ~ 1.215J 
|15j . which is almost completely accounted for by the missing energy bonds along the 
legs (0.350 J) and across the rungs (0.455 J). Then, the method developed here is a 
good approximation and can be generalized for ladders containing a low concentration 
of impurities. Below we give a simple application of this approach. Our example is done 
for classical spin systems because the calculations are almost direct in this case. Besides, 
there are also manganese halide compounds [16J that are quasi-one-dimensional and two- 
dimensional antiferromagnets. Furthermore, these Mn(II) compounds have spin 5/2 so 
they are also nearly classical and therefore, potential systems to test our results. 

At zero temperature, the Hamiltonian of spin ladders possesses, in the classical limit 
(S — > oo), a minimum given by the antiferromagnetic vacuum solution 4>q(x) = 4> z z = z, 
where z is an unit vector in the vertical direction. This solution breaks the 0(3) invari- 
ance of the model down to the subgroup 0(2) of rotations around the z-axis. Conse- 
quently there should appear two Goldstone modes, which are nothing but spin waves, 
associated with <j> x and 4> y . Then a natural first step concerning the impurity systems 
is to study the interactions between spin waves and holes. In a quantum spin system, 
the corresponding problem would be the interactions between triplons [T7] (which are 
a triplet of well defined spin-1 magnons) and holes. However, it will be considered in a 
future work. 

Expanding <f)(x,t) around the vacuum solution cj)(x,t) = <fio( x ) + v( x ^) an d min- 
imizing Hamiltonian (jUJ), one obtains, in the linearized approximation, the following 
scattering equation d^.ff(x,t) — [l/v 2 sk (x — xo)]dfff(x,t) = Uk(x)ff(x,t), where the scat- 
tering potential is Uk(x) = — {d x ln[v S) k(x — Xo)/gk(x — xo)]}d x . However, the function 
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v s ,k{x — xo)/gk{x — xq) is constant practically through all space (see Eqs,(fl~2])) while 
varies only when entering the impurity regions. Therefore, Uk{x) is zero in al- 
most all space and we have magnon solutions for the field equation in the three regions: 
(x < xq — ao), (xo — ao < x < xq + ao) and (x > xo + ao)- The form of Uk(x) is not 
explicitly known but its effects can be envisaged using the following simple analysis: if 
a magnon (for instance, coming from the left) hits the zone of influence of the potential 
Uk (x) (or the zone of the impurity) at xo — ao > then its velocity will be changed from v s 
to Vk, s = (En^fc Ja/ c ^b\Y^ 2 ano - &fter leaving behind this region at xq + ao, it will 
be changed again to v s . Consequently, supposing a plane wave coming from — oo, and 
assuming that the lowest order effect of Uk(x) is to cause elastic scattering centers for 
magnons, the solution at +00 can be approximated by tJq ex.p[i(qx—u> q t+5 nt k(q)/2)] with 
frequency co q = qv s , which of course, is precisely the dispersion relation for magnons in 
the absence of holes. The function 5 nt k(q) may be regarded as a phase-shift (for magnons 
in a ladder with n legs) which depends on the particular position (leg A;) of the hole in 
the spin ladder for a determined wave number q. Indeed, after passing the defect, the 
wave is shifted from the original one due to the different velocity Vf. s acquired in the 
hole region. This phase-shift can be easily estimated considering the difference of paths 
v s At and u^Ai, where the interval At = 2ao/ffc jS is the time necessary for the wave to 
leave behind the region of the hole (xo — ao < x < xo + ao). In the lowest order, it is 
given by 

<W(<?) = -4ga f — - 1 J = -4ga 

\Vk,s ) 

For the case of a spin ladder with two legs and N — > 00, the phase-shift does not 
depend on the position of the vacancy, which can be put at the left (k = 1) or right 
(k = 2) leg. For J\ = J2 and J a ,b = J, it is given by 

82M = <Mg) = - 4 W(i + J/2J) 1/2 - 1]. (15) 
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1/2 



^2b,c^b,c 



-1\ 1/2 
b,c \ 



(14) 



Spin ladders with three or more legs have two different situations. For example, for the 
case with three legs, if the vacancy is placed at the first or last leg, we have: 



£3,1(9) = £3,3(9) = -4<?a 



(1 + m/A) 



_(l + i?/2)(l + i?/12)_ 



1/2 



1 



(16) 



where R = J/J. By the other hand, if the hole is placed at the central leg, the phase 
shift is 



£3,2(9) = -4<?a 



(1 + 3g/4) 1 1/2 _ ' 
(l + i?/4)(l + i2/12). 



(17) 



Finally we discuss the magnon density of states in the impurity spin ladders. To do 
this, we consider a large system of size Nao and impose periodic boundary conditions 
on the continuum (magnons) states ff q (x). This periodicity, together with ff q (x) x ->± O0 s=s 
ffoex.p[i(qx ± 8 n ^(q)/2)\, gives the following condition for the allowed wave vectors: 
Na$q m + S n< k(qm) = 2ixm (m = 0, ±1, ±2, ...). Clearly, the magnon density of states is 
changed by the presence of a hole as follows 



dm Na , 1 dS n ^(q) 
en ' kiq) -^q- = ^r + ^^q—' 



(18) 



and so, the change in the density of states /S.Q n ^(q) = g n ,k(q) — Nao/2ir is given by 

1 dS n! k(q) 



(19) 



7r dq 

Above we have multiplied by a factor of 2 to take into account the two Goldstone modes. 
For a classical two-leg spin ladder one has 

-4a , 



A02,i(g) = A^ 2 ,2(9) = 



7T 



■[(l + i?/2)V2_i]. 



(20) 



For three-leg ladder 



Ag3,i(g) = Ap 3i3 (g) = 



-4o 



TT 



(1 + 312/4) 



(l + i?/2)(l + i2/12) 



1/2 



- 1 



(21) 



and 

, , -4a f T (1 + ml A) 1 1/2 1 
A ^> = V([ (l + W + fl/12) J "T P2) 

In Fig. (2) we plot -D = lS.Q n ^/ao as a function of i? = J/ J for n = 2 and n = 3. 
In general, the density of states does not depend on the wave-vector q (since the phase- 
shift has a linear dependence on q). For all legs, the limit R — ► implies D — > 0. As 
expected, it means that the presence of J is very important for the phase-shifts as well 
as for the change in the density of states. In practice, it avoids a broken lattice and 
leads to Vk, s / 0. For a three-leg ladder with the impurity placed at leg 1 or 3, Ags,i is 
very small for an appreciable range of R and becomes positive for R > 4. It increases 
considerably for large values of R and in the limit R — > oo, Ag3 i — > 1.27. Such a 
general behavior is also expected for ladders with more than three legs with impurities 
placed at the external legs (of course, it may have important qualitative changes as, for 
example, the function is positive for R < R c and then becomes positive for R > R c ). If 
the vacancy is located at the central leg, the behavior of the three-leg ladder is similar 
to the previous case but the change in the density of states is much larger. In this case, 
A^3 j 2(g)/ao becomes positive only for R > 20 (not shown in Fig.2). In addition , like 
A^3 j i(g)/ao, A^3 j 2(o , )/ao — * 1-27 as R — > oo. On the other hand, for a two-leg ladder, 
Af?2,i(<7)/oo is always negative and its modulus increases monotonically as R increases. 
The magnon density of states suffers a very expressive change for large values of J. 
Therefore, in this circumstance, there is a clear distinction in the magnon spectrum 
when a hole in the spin ladder is present or absent. 

Ladders with a higher number of legs might be approached as well, using for ex- 
ample the Mathematica program. In particular it is interesting to study how the ratio 
w = v s /vk :S changes as we increase the number n of legs and approach a two-dimensional 
lattice. We have checked this for the isotropic case J = J by putting the hole in the 
center leg of the ladder, getting the following results for n = 11,31,51,71,91 respec- 
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tively: w = 1.03358, 1.01029, 1.00605, 1.00428, 1.00331. One can immediately see that w 
approaches zero as n goes to infinity. However, the method used here for the phase-shifts 
is not convenient for two-dimensional (2d) systems because one has to deal with cylin- 
drical waves, which contains an infinite number of angular momentum channels. A more 
adequate method for the 2d case in the continuum approximation is given in Ref.|13j. 
Thus, we cannot conclude that the influence of a single static hole on the magnon phase 
shift and density of states becomes negligible in the two-dimensional system. 

We also expect that a similar behavior may happen to triplons in spin-1/2 two- 
leg ladders. Basic differences must appear due to the existence of a gap for these last 
excitations. Of course, the change in the spin wave density of states has a deep influence 
on the static and dynamical properties of ladders containing a low concentration of holes 
as it causes a change in the magnon free energy |18j . 

In summary, we have proposed an approach based on the nonlinear a model to study 
static holes in spin ladder systems with an arbitrary number of legs. As an immediate 
example of application of this method, we have studied the magnon-hole interactions in 
classical spin ladders and calculated the magnon phase-shifts and the consequent change 
in the magnon density of states. Such results may also give some insights about the 
behavior of linearized oscillatory excitations around holes in low-dimensional quantum 
spin materials. However, the calculations were performed using a long wavelength theory 
and we have to remark that in experimental situations involving small scales (of order of 
the lattice size) the application of the nonlinear a model should be viewed with caution. 
Future study considering systems containing multiple impurities (which, for adjacent 
holes, have larger sizes) and applications to quantum ladders are in progress. As in the 
two-dimensional case [HI 03]) our results may also be useful to study possible topological 
excitations interacting with holes in spin ladders. 
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Figure Captions 

Figure 1. Configurations of spin ladders along the z-direction with two and three legs 
containing a hole. For the case of three-leg ladders, the second and third configurations 
are equivalents, but the fourth leads to different results. 

Figure 2. Change in the density of states D = Ag n ^/ao as a function of the rung 
coupling R = J I J for spin ladders with two-leg (solid line) and for the two possibilities 
for three-leg: vacancy placed at the central leg (dashed line) and at leg 1 or 3 (dotted 
line). 
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